Exact solutions for mass-dependent irreversible aggregations 
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We consider the mass-dependent aggregation process (k + l)X — > X, given a fixed number of unit 
mass particles in the initial state. One cluster is chosen proportional to its mass and is merged into 
one either with fc-neighbors in one dimension, or - in the well-mixed case - with k other clusters 
picked randomly. We find the same combinatorial exact solutions for the probability to find any 
given configuration of particles on a ring or line, and in the well-mixed case. The mass distribution 
of a single cluster exhibits scaling laws and the finite size scaling form is given. The relation to the 
classical sum kernel of irreversible aggregation is discussed. 
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Recently the theory of irreversible aggregation was 
revisited in view of renormalization of complex net- 
works [![. In Q, a simple mapping between random 
sequential renormalization (RSR) [2|, |3j and irreversible 
aggregation Q was pointed out, where a local random 
renormalization step to produce a new 'super-node' in 
complex networks corresponds to one aggregation event 
of 'molecules'. Exact combinatorial analyses, both in 
one dimension (without diffusion) and in the well-mixed 
case, gave the same scaling law of cluster mass distribu- 
tion. The corresponding exponent only depends on k, 
the number of interacting neighbors This RSR pro- 
cedure corresponds to the 'constant' kernel of irreversible 
aggregation (one of three well-known "classical" kernels- 
constant, sum, and product kernels [4]). In this study, 
we show the relation between mass-dependent RSR and 
irreversible aggregation with the sum kernel. Applying 
the same combinatorial technique of [l[ , we find the exact 
solutions for mass-dependent irreversible aggregation as 
well. 

Here we consider models governed by the reaction 
(k + 1)X — > X, where a cluster is picked randomly, in 
proportion to its mass, after which it coalesces with k 
other clusters. In the case of one-dimensional models 
these are fc-neighbors, while they are k other clusters cho- 
sen randomly in the case of well-mixed systems. In both 
cases, the other clusters are chosen independent of their 
masses. The mass of the newly formed cluster is the sum 
of the (k + 1) masses. For one dimensional models, both 
a ring with periodic boundary condition and a line with 
open boundary condition are considered. Reactions are 
allowed only if there is a sufficient number, k, of available 
clusters. 

First, let us consider the model defined on a 'ring'. Ini- 
tially, Nq particles of unit mass (to = 1) are placed on 
a ring like beads (see Fig. []}. Each particle is labelled 
by i <E [l,...,No\. At each time, one cluster is picked 
in proportion to its mass, and is subsequently merged 
with its k right neighbors into one big cluster having a 
mass equal to the sum of the (k + 1) masses. Cluster 
masses are therefore restricted to m = 1 (mod k). This 
can be written as m — 1 = ks, where s is the number of 



aggregation events needed to make a cluster of mass to. 
Similarly the number of clusters at any time, t, is given 
by N = No — kt, where time t is denoted by positive inte- 
gers representing the total number of aggregation events. 
We do not allow two events to happen simultaneously 
in this study. Otherwise they can happen either at reg- 
ular intervals, intermittently, or according to a Poisson 
process. 

To find the probability that any of the N clusters 
picked at random has mass to resulting from s aggrega- 
tion events, we follow an approach similar to the one in- 
troduced in Ref. [1]. The crucial observation that makes 
the analysis simple is that picking clusters according to 
their mass is equivalent to picking sites with uniform 
probability, since a cluster of mass to occupies m sites. 
Let i be any site {e.g. i — 1), and let 7t^°(to) be the 
probability that a cluster of mass m starts at this site 
and occupies the sites (i, i + 1, . . . ,i + m— 1). The prob- 
ability that any of the N clusters picked at random has 




FIG. 1: (Color online) Illustration of aggregation on a ring 
with k — 2, No = 24, and N = 6. The tree in color corre- 
sponds to a cluster of mass m = 7. It has seven leaves (red) 
and three internal nodes (blue). The numbers beside internal 
nodes correspond to the time when coalescence occurs. 
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mass m after t events is then 
and 



(m), 



N, 



(m) = 
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^total 
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where "cluster is the number of possible histories of aggre- 
gation events i 2 , . . . , i s ) leading to a cluster of mass 
m, "rest is the number of possible ways to form the other 
(JV— 1) clusters, and "total is the total number of histories 
for t merging events. The binomial coefficient (*) corre- 
sponds to the number of choices associated with different 
time orderings for the s events in the cluster of mass m 
and the (t — s) events in the rest of the clusters. 

The total number of all histories involving t events is 
simply 



"total = N Q . 



(3) 



This is to be contrasted to the number of histories n^ ,JV °' 
that lead to the first cluster starting at i = 1 and the N- 
th ending at Nq. A somewhat more involved argument 
gives 



(4) 



The number of histories leading to a single cluster of mass 
m covering the sites of interval [1 , m] is thus 



_ [l.m] _ s _l 
"cluster — "i — m i 

while 

"rest = 4'5r ml ={N-l)x (N - m)'-*- 1 
Combining Eqs. (jT|)— (J6j) , we finally obtain 

N - 1 ft\ m s - x (iVo - my- 3 - 1 



PN°( m ) = 



N 



N' 
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For this mass-dependent aggregation process, we can 
also work out the joint probability distributions for 
masses of adjacent clusters. We denote by p%°(mi,m2) 
the probability to find a cluster of mass m\ followed im- 
mediately to the right by a cluster of mass m 2 . This is 
non-zero only if m\ = ks\ + 1 and m,2 = ks 2 + 1, where 
s a is the number of aggregation events needed to form a 
cluster of mass m a . By the previous arguments, we get 



Px°(rnx,m 2 ) 



N - 2 
N 



t 



where s = t — X)s=i s ana - TO o = -^o — Ss=i m P- ^ 
is interesting to note that this joint probability properly 
holds the following relation, 



For any 1 < a < N — 1, the joint probability distribu- 
tion for a consecutive adjacent clusters is given by 



Pm (mi. —,m a ) 



N-a T[t, {s},a + 1] n^=o T77 - 
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where we used the multinomial coefficient 



(8) 



T[t,{s},a+1] = 



Sq,...,S c 



In particular, this can be done for the joint distribution 
for all N masses by setting a = N — 1. The resulting 
expression is always invariant under any permutations 
of N numbers (mi, ...,mjv}, as was the case with mass- 
independent aggregation Hence the iV-cluster prob- 
ability is independent of the spatial ordering of the clus- 
ters. Therefore, there are no spatial correlations, even 
though there are obvious correlations between the masses 
at any given time. For this reason (and as verified in de- 
tail using Eq. (QJ instead of Eq. (|3J|), the joint probability 
for N masses on a line, i.e., a one-dimensional system 
with open boundaries, is also given as Eq. ©, showing 
that the two models lead to the same statistics for any 
o. 

The absence of spatial correlations indicates that the 
same dynamics might also result from the well-mixed 
case. To check this, we now start with a bucket con- 
taining Nq balls, each of unit mass. An event consists 
of first picking one ball with probability proportional to 
its mass and then choosing k balls out of the bucket, in- 
dependent of their masses. The balls are merged and a 
new ball, having a mass equal the sum of the masses of 
its (k + 1) constituents is returned to the bucket. This 
process repeats until N clusters remain. 

The single cluster mass distribution for the well-mixed 
model can be obtained using the same strategy as before. 
Since events now correspond to choosing one ball with a 
mass- weighted probability, and k balls out of (No — kt—1) 
balls randomly, we have a i-power of Nq and a product 
of binomial coefficients, 



. r ,'JVb - A A T [N + k - 1 
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The expressions for n c i ustG r and n 1GSt are analogously 
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The number of time orderings is exactly the same as 
before, but the first factor Nq/N in Eq. ([I} has to be 



3 



z° z 
cl 



10 u 



10" 1 



10" 



10- 




□ N=2 
o N=6 
a N=10 
v N=14 
— exact solutions 



10 u 



10' 
m 



1(T 



FIG. 2: (Color online) Cluster size distributions after t — 49 
(N = 2) events for k — 2, for different values of AT averaged 
over 10 6 realizations compared to exact results. The power- 
law slope of small m is —3/2 independent of k. The large size 
behavior changes from an increasing power law to a decreas- 
ing one around N ~ V No. The inset shows the symmetric 
distribution for N = 2. 



replaced by (^j 5 ) . Putting all these considerations to- 
gether, many cancellations take place, leading exactly to 
Eq. (0) again. This argument can be similarly extended 
to get the full iV-particle distribution function, obtaining 
exactly the same result as Eq. (jHJ for any k and a. 

Let us look at the characteristics of the solutions. Even 
though the composition principle is the same as that of 
mass-independent aggregation in Ref. [H, the final so- 
lution and the characteristics are quite different. First, 
cluster size distributions at several different times are 
shown in Fig. [2] for N a = 100 and k — 2. The symbols 
indicate the numerical simulation results over 10 real- 
izations and the solid lines are the exact solutions of Eq. 
([7]). The tail region corresponding to large cluster sizes 
changes from a fast exponential decay to an increasing 
power law as the merging process approaches termina- 
tion. The turning point is around N ~ y/T^o- When 
N = 2, since the sum of the two cluster sizes is always No, 
the distribution p^° (to) is symmetric under the exchange 
to O N — to for any k. The symmetric distribution for 
N = 2 is shown in the inset of the Fig. [2l 

When No — > oo, asymptotic power laws can be deter- 
mined using Stirling's formula. If N is fixed and both m 
and (No — to) —¥ oo, one obtains the scaling form 



p%o(m) ~ N 



N 2 f 




N o, 



(12) 



For small masses, this gives a decreasing power law, with 
exponent —3/2, independent of k. Interestingly, this 
is very different from mass-independent aggregation, for 



which the analogous exponent depends on k and is equal 
to — 1 + 1/fc [lj. The exponent —3/2 is the same as that 
for the aggregation with the sum kernel of the irreversible 
aggregation obtained in Ref. Indeed, the rate equa- 
tion for the current aggregation model in mean-field the- 
ory is the same as for the sum kernel 4, 5]. For the k = 1 
case, the rate equation is simply 

m oo 

Ap m = ^2 m'Pm'Pm-m' ~ Pm ^ (m + m')p m > 

m' — 1 m' — l 

m 

= — ^ P m >Pm-m' ~ (m + fh)p m , (13) 



m' — l 

where p m denotes p^°(m) to make the equation more 
concise, and where to means the mean cluster size. Equa- 
tion (|13[) is the same as the rate equation for the sum 
kernel in [1, @. The behavior for large to is different, 
however, and is not described by mean field theory. 

According to Eq. (fT2|) . N Q 2 p^° (to) should be a function 
of m/No only for fixed 



y 



N 

kN~o 



(14) 



and for N <C A^o, No — to. The resulting data collapse is 
shown in Fig. [31 where we also factored out a power of 
m/No to make the curves less steep. Notice that N and 
No are related by N = N (mod k), which implies that 
the values of y used in this plot are not strictly constant 
but deviate slightly from their nominal values for small 
No, which causes the deviation from a perfect collapse 
for y = 0.2. Even for different values of k, this scaling 
function works, as can be seen in Fig. U where the scaling 
collapses for three cases, k = 1,2, and 4, are shown. Sur- 
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FIG. 3: (Color online) Finite size scaling collapses for different 
y and fixed k — 2. Collapse lines for y = 0.2 and y = 1.0 are 
shifted up and down to make them distinguishable from other 
collapse lines. 
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FIG. 4: (Color online) Finite size scaling collapses for different 
k. In order to check the scaling collapses for different k, kNo 
is fixed at 400 and TV = 4, 16, and 24, which correspond to 
y = 0.04, 0.64, and 1.44. Collapse lines for y = 0.04 and y = 
1.44 are shifted up and down to make them distinguishable 
from each other. 



prisingly, this means that the process of choosing a clus- 
ter proportional to its mass in conjunction with choosing 
two clusters at random for k = 2 is asymptotically the 
same as repeating the merging process for k — 1 twice in 
the sense of the scaled mass. 

The probability p^ a (m) satisfies the following recur- 
sion relation 



N ~N+1 , 

p%> +k (m) = A J2' ^?(4hiW (15) 

m' =m-\-k 



with 



.4 



N N(k + l) 
(N -N)(N + ky 



(16) 



where the prime on the summation symbol indicates 



that 



must increase in steps of k. Interestingly 



this quadratic recursion relation corresponds to the 
time-reversed process of aggregation, i.e., fragmenta- 
tion. As with the quadratic recursion relation of mass- 
independent aggregation [l[, the mass distribution at 



N + k is given by the product of the mass distribution 
at N describing the relative probabilities with which the 
cluster fragments, given by p™ +1 (m), and the total frag- 
mentation probability. The latter was just oc (m! — 1) 
in the mass- independent case [l| , while now it is propor- 
tional to m 7 1 ■ Equation (|15p follows then by considering 



how fragmentation leading to a cluster with mass m goes 
through an intermediary with mass m' . 

We also examined numerically the aggregation pro- 
cesses where the clusters were chosen with probabilities 
proportional to higher powers of their mass, in particular 
oc m? and oc m 3 , i.e., the square of a cluster's mass and 
the cubic of a cluster's mass. The asymptotic power law 
exponents are roughly about —5/2 and —7/2 for the m 2 - 
dependence and m 3 -dependence respectively. However, 
exact solutions for these cases have not yet been found. 

In summary, we derived the exact solutions for the 
probabilities to find any configuration after a fixed num- 
ber of aggregation events in the models where a cluster 
picked with probability proportional to its mass aggre- 
gates with k other particles. More specifically, we studied 
three versions of this process (particles on a ring join- 
ing with nearest neighbors, particles on an open-ended 
line, and the well- mixed case), and found exactly the 
same solutions using combinatorial counting. We at- 
tribute this to the absence of spatial correlations, al- 
though they are a priori not excluded. Differently from 
the mass- independent random sequential renormalization 
(RSR), which shows fc-dependent exponents in scaling 
laws for small masses, the cluster size distribution fol- 
lows a power-law with exponent —3/2 independent of k, 
which is the same with that of the sum kernel for irre- 
versible aggregation. Finally, the aggregation process is 
also related to a time-reversed fragmentation process, the 
characteristics of which are briefly discussed. 

Mass dependent RSR and the related aggregation pro- 
cess was also considered in two dimensions Q, where a 
runaway giant cluster exists after few steps and takes all 
merging action. The behavior is very similar to the gela- 
tion in the aggregation process with the product kernel, 
but aggregation events in two dimensional RSR involve 
fluctuating numbers of neighbors, differently from the ag- 
gregation process considered in the present paper. 
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